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The agreement of hydrodynamic predictions of differential elliptic flow and radial flow patterns 
with Au+Au data at ^snn = 200 GeV is one of the main lines of evidence suggesting the nearly 
perfect fluid properties of the strongly coupled Quark Gluon Plasma, sQGP, produced at RHIC. We 
study the sensitivity of this conclusion to different hydrodynamic assumptions on hadro-chemical 
and thermal freezeout after the sQGP hadronizes. We show that if chemical freezeout occurs at 
the hadronization time, as required to reproduce the observed hadron yields, then, surprisingly, the 
differential elliptic flow, V2(pr), for pions continues to increase with proper time in the late hadronic 
phase until thermal freezeout and leads to a discrepancy with the V2(pr) data. In contrast, if both 
hadro-chemical and thermal equilibrium are maintained past the hadronization point, then the 
mean transverse momentum per pion increases in a way that accidentally preserves V2(pr) from the 
sQGP phase in agreement with the data, but at the cost of the agreement with the observed hadronic 
yields. In order that all the data on (1) hadronic ratios, (2) radial flow, as well as (3) differential 
elliptic flow be reproduced, the sQGP core must expand with a minimal viscosity, r\ m that is 
however even greater than the viscosity, t\h ~ T/o~h, of its hadronic corona. However, because of 
the large entropy density difference of the two phases of QCD matter, the larger viscosity in the 
sQGP phase leads to nearly perfect fluid flow in that phase while the smaller entropy density of 
the hadronic corona strongly hinders the applicability of Euler hydrodynamics in that phase. The 
"perfect fluid" property of the sQGP is thus not due to a sudden reduction of the viscosity at the 
critical temperature T c , but to a sudden increase of the entropy density of QCD matter and is 
therefore an important signature of deconfinement. 



PACS numbers: 24.85.+p,25.75.-q, 24.10.Nz 

I. INTRODUCTION 

One of the most intriguing experimental findings at the 
Relativistic Heavy Ion Collider (RHIC) in Brookhaven 
National Laboratory (BNL) is the large magnitude of 
the elliptic flow parameter vi [1-3] in comparison with 
the smaller values observed at lower collision energies 
(for results at Super Proton Synchrotron (SPS) energies, 
see Refs. [4-6]). The magnitude of v-i and in particular 
its transverse momentum px and mass m dependences 
at RHIC were found to be close to predictions based on 
ideal, non-dissipativc hydrodynamics simulations around 
midrapidity (| r\ | < 1), in the low transverse momentum 
region (pr^S 1 GeV/c), and up to semicentral collisions 
(b< 5 fm) [7, 8]. This result has led to the recent BNL 
announcement [9] about the discovery of the near perfect 
fluidity of the strongly coupled/interacting quark gluon 
plasma (sQGP) [10-12] produced in ultra-relativistic nu- 
clear reactions at RHIC. 

Until RHIC data, "perfect fluidity" was never ob- 
served nor expected to apply theoretically in high en- 
ergy hadronic or nuclear reactions due to nonvanishing 
viscous dissipation [13]. Especially, since the discovery of 
asymptotic freedom in QCD, the prevailing paradigm has 
been the expectation of large viscosities in a weakly cou- 
pled/interacting QGP (wQGP) at very high densities. In 
addition, it is well established [14] that the hadronic reso- 
nance gas phase of QCD matter is highly dissipative. The 



discovery of elliptic flow at RHIC consistent with nearly 
perfect fluidity is therefore an experimental and theoreti- 
cal surprise. Hence a new name, sQGP, has been adopted 
to characterize the observed strong coupling properties 
of the QGP near the critical temperature T c ~ 160-170 
McV that keep viscous effect to a minimum at RHIC. 

In this paper, we present the case for the following 
physical interpretation of RHIC data based on current 
hydrodynamic analyses: (1) the high density core part 
of matter produced in relativistic heavy ion collisions, 
i.e. the sQGP, must expand as a nearly perfect fluid 
despite of its higher viscosity, (2) the perfect fluidity of 
the sQGP core is a consequence from a large jump of 
the entropy density at the critical temperature, T c , i.e. 
deconfinement, and not from some anomalous reduction 
of its viscosity, (3) viscous effects on its hadronic corona 
are necessarily large despite its smaller viscosity, and (4) 
ideal inviscid hydrodynamics should not be applied to the 
hadronic corona which requires a noncquilibrium trans- 
port description. 

In Sec. II, we discuss why we expect a surprising mono- 
tonic increase of the viscosity of QCD matter through the 
critical temperature and emphasize the important role 
played by the rapidly varying viscosity to entropy ratio 
in connection with perfect fluidity of the sQGP phase. In 
Sec. Ill, we discuss different assumptions for the hadronic 
matter in the hydrodynamic models to clarify what are 
open issues in the current hydrodynamic approaches. In 
Sec. IV, the time evolution of the transverse energy per 
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particle is discussed. The mean transverse energy is 
found to be the key to distinguish the model assumptions 
in the hadron phase. Results from the hydrodynamic 
simulations are reviewed in Sec. V. We will show how 
the perfect fluid description for the hadronic matter in 
chemical equilibrium in the conventional hydrodynamic 
simulations leads to accidental reproduction of px spec- 
trum and V2(pt)- In order to understand analytically the 
role of chemical freezcout on the transverse dynamics, we 
employ a blast wave model and give a dynamical mean- 
ing to this model in Sec. VI. Finally, summary of this 
study and an outlook are presented in Sec. VII. 



II. VISCOSITY AND ENTROPY IN QCD 

Weak coupling pcrturbative QCD (pQCD) estimates 
[15-17] of the viscosity of a wQGP were based on basic 
kinetic theory relations 
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iTs SB {T) » 3Tn SB (T) « K SB T 4 is the energy density, 
pressure, entropy density, and number density of an ideal 
Stefan-Boltzmann (SB) gas of quarks and gluon charac- 
terized by the constant A' SB = f^[2(JV c 2 - 1) + |l2iV/] ~ 
12-15 for N c = 3, Nf = 2-3. The key microscopic dy- 
namical quantity in Eq. (1) is the transport mean free 
path Atr = l/("sB0tr) which is controlled in pQCD by 
the Debye screened transport cross section [16, 18] 



Otr 



= / ddei s'm 2 cm 



"(1 + 2) 



(2z + 1) In 1 + 



, (2) 



where z = /i 2 /s and s f=a 17T 2 is the mean partonic 
Mandelstam variable. Perturbatively, the screening mass 
squared varies as /z 2 ~ Aita s T 2 . For numerical estimates 
we take a s (T) w 4tt/[18 ln(4T/T c )], so that a s (T c ) ~ 0.5. 
In the range T c < T < 5T C (0.5 > a s (T) > 0.23), the 
perturbative transport cross section cr t r < 2 mb remains 
much smaller than typical hadronic cross sections <jh ~ 
10-20 mb [16, 19, 20]. 

An important dimcnsionless measure of how imperfect 
or dissipative a fluid may be given by the ratio of vis- 
cosity to entropy density, rj/s [21]. This is most easily 
seen via the Navier-Stokes equation in (l+l)-dimensional 
boost invariant hydrodynamics [15, 16, 19]. In the perfect 
(Eulcr) fluid limit, the proper energy density decreases 
with proper time, r, due to longitudinal expansion dV = 
■KR 2 dr and PdV work via de/dr = -{e + P)/r = -sT/t 
with a solution T = ^(tq/t) 1 / 3 for massless particles 
[22]. However, shear and bulk viscosity reduce the abil- 
ity of the system to perform useful work by adding a 



term (477/3 + C)/ 7 " 2 to the right hand side. Neglecting 
bulk viscosity, that vanishes in equilibrium for massless 
partons, de/dr = -(sT/r)[l - 4(?//s)/3Tt]. This shows 
that for the earliest times consistent with the uncertainty 
principle [16], r ~ 1/3T, the cooling of the plasma due 
to both expansion and work is canceled if rj/s > 1/4. 
The ability of the system to convert internal energy into 
collective flow is thus severely impaired at early times 
if rj/s and £/s are not very small. In fact, in order to 
reproduce the observed elliptic flow at RHIC, numerical 
solutions to covariant parton transport equation [18] and 
blast wave analysis with viscous corrections [23] showed 
that rj/s had to be less than about 0.2 during the first 3 
fm/c. 

The viscosity to entropy ratio in the weakly coupled 
QCD plasma on the other hand is 
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This ratio is not small (n/s = 0.35,0.48,0.58,0.66 for 
T/T c = 1,1.5,2.0,2.5) indicating that the wQGP is ex- 
pected to be a rather "poor fluid" with large dissipative 
corrections. 

The analytic dependence on a s in Eq. (3) reproduces 
well the approximate temperature dependence implied 
by Eq. (2) if we assume the perturbative variation of 
the screening scale ji. Lattice QCD data [24] indicate 
that /x ~ (2.0-2.5)T is not far from the perturbative es- 
timate extrapolated into the physical g > 1 region and 
that ot s (T) < 0.5 above T c is also reasonable. However, 
the underlying simple gas kinetic approximation for the 
viscosity is only rigorously valid in the g <C 1 region. 

Formally, by increasing a s > 0.5, it would seem that 
the right hand side of Eq. (3) could be made to be as small 
as we like if we ignore the ln(l/a s ) singularity. However, 
by the Heisenberg uncertainty principle, the transport 
mean free path cannot be localized to less than Ax ~ 
l/(p) ~ 1/3T. This leads to a quantum kinetic lower 
bound on the viscosity for ultrarelativistic gases [16]: 
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with an undetermined multiplicative factor on the order 
of unity. 

A special quantum field theoretic determination of a 
viscosity lower bound was found recently for infinitely 
coupled supersymmetric Yang-Mills (SYM) gauge theory 
using the Anti de-Sitter Space/Conformal Field Theory 
(AdS/CFT) duality conjecture [25]: 
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This bound is obtained in the dual N c = 00 and g 2 N c = 
00 limits of the special Af = 4 conformal SYM [26]. It 
is interesting to note that this analytic SYM bound is 
close to the simple kinetic theory uncertainty principle 
bound in Eq. (4). It has been conjectured [25] that 1/4-7T 
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in Eq. (5) is the universal minimal viscosity to entropy 
ratio even for QCD. In that case, the viscosity of the 
sQGP could be up to a factor of ~ 1/2-k smaller than 
of a wQGP. It is then tempting to conclude that the 
sQGP must have anomalously small viscosity if perfect 
fluid behavior is observed. However, as we show below, 
the sQGP viscosity is actually very close to that of ordi- 
nary hadronic matter just below T c . 

To develop this argument further, we first digress to 
recall that the entropy density in the N c ^> l,g 2 N c ^> 1 
limits of TV = 4 SYM is given by [27] 
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where the Stcfan-Boltzmann constant for AT = 4 SYM is 
Ksym = 7I ' 2 (-^c — ~ 39.5 is about 3 times greater 
than K$b of our QCD world [26]. What is especially re- 
markable about Eq. (6) is that, at infinitely strong cou- 
pling, the entropy density is only reduced by ~ 25% from 
its non-interacting SB value. On the other hand, the vis- 
cosity in this extreme limit is reduced about an order 
of magnitude from the weak coupling value and limited 
only by the quantum (Heiscnberg uncertainty) bound on 
the effective scattering rate. Current lattice data on the 
QCD viscosity near T c [28] are with large numerical error 
bars between these weak and super strong coupling limits 
but the relatively small deviation of the lattice entropy 
density from the SB limit is consistent with Eq. (6). 

The AdS/CFT lower bound (5) together with the as- 
sumed universal 3/4 reduction of the SB entropy density 
implies that the absolute value of the sQGP viscosity at 
T r would be 
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where we used a fact that for QCD K$b ~ 12-15 is acci- 
dentally close to 47r. The monotonic increase of ?7sym(T) 
is illustrated by the dashed curve in Fig.l. 

The effective transport cross section via Eq. (1) at T c ~ 
160 MeV is in this case 
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Here r] c = T c 3 = 0.106 GeV/fm 2 at T c = 0.16 GeV. See 
Ref. [29] for an independent estimate of the transport 
cross section in the sQGP phase leading to similar a^iT) 
near T c . 

While there is no consensus yet on what physical mech- 
anisms could enforce the minimal viscosity bound in the 
sQGP [18, 30, 31], we take as empirical fact that the 
sQGP viscosity must be close (within a factor of two) 
to the minimal (uncertainty) bound, Eq. (7). Our cen- 
tral assumption is that local thermal equilibrium is main- 
tained in the sQGP core with minimal dissipative devia- 
tions and with the equation of state and hence speed of 
sound as predicted by QCD. Alternate scenarios, with ar- 
bitrary equations of state with higher speed of sound that 
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FIG. 1: Illustration of the approximately monotonic increase 
of absolute value of the shear viscosity with temperature. 
The kink shown at T c is expected to be smeared out by the 
AT C /T C ~ 0.1 width of the QCD cross-over transition. The 
solid blue curve shows rj(T < T c ) = T/gh for a HRG followed 
by the more rapid increase of the viscosity in the sQGP phase 
with ??sQGP ~ ?7sym = KsbT 3 /A-k w T z . The horizontal line 
shows that near T c , r\ « r\ c = T 3 . At high r>T c asymptotic 
freedom leads to an even more rapid growth of viscosity as 
the sQGP evolves gradually into the weakly coupled wQGP. 
In this figure, w = 1 in Eq. (10) is taken to emphasize the 
possibility that the highly viscous but nearly "perfect fluid" 
sQGP may become an ordinary "viscous fluid" already for 
T > 2T C . 



in principle could compensate the higher dissipation and 
viscosity in a wQGP will not be considered here. In this 
connection we also emphasize the importance of fixing 
sQGP initial conditions with Color Glass Condensate or 
saturating gluon distributions constrained by the global 
entropy observables [11, 32]. With fixed initial conditions 
and equation of state, the remaining degrees of dynam- 
ical freedom are reduced to the dissipation corrections 
discussed in this section for the sQGP phase and the 
dynamical constraints on its dissipative hadronic corona 
discussed in the subsequent sections. 

Note that the effective transport cross section in the 
sQGP of r just above T c is remarkably close to the 
hadron resonance gas transport cross section just below 
T c [19, 20]. However, due to the 1/T 2 scaling at T ~ 2T C , 
the effective transport cross section in the sQGP would 
already drop to ~ 3 mb while preserving the (uncertainty 
principle) lower bound Eq. (5). 

In contrast to the novel sQGP phase above T c . for 
T < T c , matter is well known to be in the confined hadron 
resonance gas (HRG) phase where the kinetic theory vis- 
cosity [16, 19] is 
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as illustrated by the solid curve below T c in Fig. 1. Be- 
cause the hadronic transport cross sections are typically 
ajj ~ 10 — 20 mb, the combination of Eqs. (7) and (9) 
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FIG. 2: Illustration of the rapid variation of the dimension- 
less ratio of the shear viscosity, Tj(T), to the entropy den- 
sity, s(T). The sharp discontinuity illustrated is not due to 
a rapid change of the transport coefficient (see Fig. 1) but to 
the rapid increase of the entropy density in QCD near T c . As 
in Fig. 1, we expect the discontinuity to be smeared into a 
rapid drop within AT C /T C ~ 0.1. Solid (dashed) blue curve 
illustrates the change of 77/ s of a HRG with c 2 H — 1/3 (1/6), 
sq/sh = 10 (3) into an approximate "perfect fluid" sQGP 
at T c . The red long dashed curve is (?7/s)sym = 1/47T. At 
T>r t asymptotic freedom gradually transforms the sQGP 
into an ordinary viscous fluid wQGP (green) , here shown for 
w= i 1. 



shows that we should not expect a large variation of the 
absolute value of the matter viscosity across T c if the 
minimal rj/s holds above T c . In Rcf. [19], Gavin found 
that for a pion gas with P-wave p and D-wave /° reso- 
nance interactions, the thermal averaged transport cross 
section and viscosity from his Fig. 3 are (r]/r] c , 07/ [mb]) ~ 
(0.66,20), (0.9,17), (1.1, 15) for T = 180, 160, 140 MeV. 
In Ref. [20], Muronga used the UrQMD resonance gas 
Monte Carlo to compute rj(T)/r] c ~ 0.75, 1.1, 1.9 for 
T = 0.14, 0.16, 0.18 GcV. In both studies [19, 20] nu- 
merical estimates are thus consistent with Eq. (9) for 
T < T c . For nonvanishing baryon density, see recent es- 
timates in Ref. [33]. 

In the sQGP phase, the minimal viscosity, Eq. (7), is 
predicted to grow cubically with T beyond T c . However, 
at T T c asymptotic freedom predicts that it would 
grow even more rapidly as the sQGP transforms grad- 
ually into a wQGP. An interpolation formula between 
these phases can be constructed as 



ri(T) 



T 3 1 (T/r c y 



T<T C 



(T/T c ) 3 [l + w(T) ln(T/T c )] 2 , T > T c 



(10) 

The extra squared logarithmic growth of the viscosity at 
T ^> T c is expected from Eq. (3). To be consistent with 
the perturbative wQGP at T 3> T c one should take 
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possible scenario may be that w ~ 1 already near T c . 
This possibility is shown in Fig. 1 by the solid curve 
above T c which would imply sQGP^ wQGP already 
above ~ (2 — 3)T C . In fact, current lattice data on the 
evolution of screening scales in the QGP phase suggest 
that hadronic scale correlations may persist only up to 
T ~ 3T C [34, 35]. A value w(T > 2T C ) ~ 1, is also not in- 
consistent with current lattice results [28]. We note that 
future measurements of elliptic flow in A+A collisions at 
LHC with ^snn = 5500 GeV will be able to test exper- 
imentally if such a precocious onset of dissipative wQGP 
dynamics occurs. 

The approximate continuity of the viscosity across T c 
indicated in Fig. 1 is to be understood to hold within a 
factor on the order of unity. What changes rapidly at T c 
is not the viscosity of QCD matter but rather its entropy 
density due to the dcconfincmcnt of the quark and gluon 
degrees of freedom. 

For a hadronic resonance gas charactered by a speed of 
sound cjj, the entropy density s(T) — sh{T /Tc) 1 ^ 11 with 
decreasing temperature decreases much more rapidly 
than does the viscosity for typical c 2 H ~ 1/6-1/3. Just 
beyond T c -possibly only up to several times T c — it is pos- 
tulated that the sQGP phase may exist with rj/s < 0.2 
but with viscosity close to ideal gas. 

Summarizing the discussion up to now, we expect that 
77 varies smoothly near T c as in Fig. 1 while the ratio 
rj/s may have a significant discontinuity due to the rapid 
onset of deconfinement as a function of T as shown in 
Fig. 2. We therefore propose the following interpolation 
formula for the temperature dependence of the rj/s ratio 
with a T independent constant w 
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We illustrate Eq. (12) in Fig. 2. Note that it is the en- 
tropy jump sq/s h ~ 3-10 that causes a drop of r]/s 
across T c . Since the HRG— ► sQGP transition with dy- 
namical quarks appears from the lattice QCD to be only 
a rapid crossover, the discontinuity is understood to be 
spread out over a temperature range AT C /T C ~ 0.1. 



III. IMPERFECT FLUIDITY OF THE 
HADRONIC CORONA 

In the last section we presented the case that the rj/s 
ratio may be small enough above T c in the sQGP core for 
perfect fluidity to arise during the critical first ~ 3 fm/c, 
while the spatial azimuthal asymmetry of the matter pro- 
duced in non-central reactions is large enough to induce 
collective elliptic flow. However, during the subsequent 
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~ 10 fm/c evolution after hadronization ol the sQGP 
core, the whole system evolves as HRG corona. In the 
HRG, rj/s is too high for local equilibrium to be main- 
tained due to its small entropy density compared with 
sQGP. Nevertheless, the data on V2(pt) seem to agree 
very well with some hydrodynamic predictions based on 
the assumption that local equilibrium is maintained until 
thermal freeze-out. However, various assumptions about 
the hadro-chemical evolution are known to lead to sig- 
nificantly different predictions for the differential elliptic 
flow. In this section we focus on the question of the va- 
lidity of the application of hydrodynamics to analyze the 
entire evolution in A+A at RHIC. 

The key problem that we now address is the role of fi- 
nal state hadronic interactions in possibly modifying con- 
clusions inferred about the prefect fluidity of the sQGP 
core. In order that the sQGP elliptic flow signature of 
perfect fluidity survives during the evolution through the 
extended hadronic "corona" we must study how longitu- 
dinal flow, transverse radial flow, as well as the elliptic 
deformation of the transverse flow may evolve in hadronic 
matter. 

Several puzzling features suggest the importance of 
investigating more closely the distortions that can be 
caused by final state hadronic interactions involving 
hadro-chemical and thermal freezeout. In ideal hydro- 
dynamics it is well known that while the central rapidity 
region is well reproduced by hydrodynamics, this is not 
the case in forward/backward rapidity regions [8]. Hy- 
drodynamics also strongly overestimates vi at energies 
below ^snn =200 GeV as well as in the most peripheral 
collisions where initially a larger fraction of the trans- 
verse elliptic interaction region starts out in the hadronic 
phase. All these data point to the fact that the dynamics 
in the hadronic corona is not at all ideal. 

Another important issue in ideal hydrodynamics as 
well as in other dynamical models is the so-called HBT 
puzzle [36] . In spite of the apparent success of hydrody- 
namic description for elliptic flow, hydrodynamics fails to 
reproduce the experimental data of the HBT radii [37- 
39]. The best current description of hadron freezeout 
consistent with the HBT puzzle involves an assumption 
of a highly dissipative resonance gas dynamics and trans- 
port [40]. 

As compiled in Fig. 20 in Ref. [41], some hydrodynamic 
results reproduce neither «2(pt) nor pt spectrum. This 
immediately raises the following two questions: 

(Ql) Are hydrodynamic results consistent with each 
other at RHIC energies? What assumptions lead to the 
differences among hydrodynamic predictions? 

(Q2) How robust is the statement that hydrodynamic 
description at RHIC data is good at low pt? 
The differences arise from the treatment of hadron phase 
dynamics. The treatment of the sQGP phase is essen- 
tially the same in all approaches so far: Parton chemi- 
cal equilibrium and inviscid hydrodynamics are assumed 
in the sQGP phase. There are, however, three generic 
classes of approaches to the evolution of the hadronic 



corona in the literature. 

Chemical equilibrium model (CE). Most of the hydro- 
dynamic calculations so far are based on the assump- 
tion that the hadron phase is a perfect fluid in both 
hadro-chemical and thermal equilibrium. With this as- 
sumption, the centrality, transverse momentum, and/or 
(pseudo)rapidity dependences of elliptic flow are studied 
[7, 8]. However, it is known that the yields of heavy 
hadrons (essentially all hadrons except for pions) are 
smaller in CE than data since the numbers of particle arc 
counted on the hypersurface at thermal freezeout within 
this approach. At rclativistic collisional energies, ther- 
mal freezeout temperature T th is smaller than chemical 
freezeout temperature T ch [42, 43]. So the numbers of 
heavy particles are suppressed due to the Boltzmann fac- 
tor and, eventually, lead to discrepancy from the data. 
CE therefore fails to account for the observed particle 
abundance systematics [44]. 

Partial chemical equilibrium model (PCE). In 
Refs. [45-48], chemical freezeout being separated from 
thermal freezeout is taken into account in the hydrody- 
namic simulations toward simultaneous reproduction of 
particle ratios and particle spectra. Below T ch , one intro- 
duces chemical potential for each hadron associated with 
the conserved number. Inelastic processes only through 
strong interactions are supposed to be equilibrated in 
the hadron phase, e.g. ji p = 2p^, /iA = Pn + Ptt, and 
so on. Note that the conserved pion number within 
this approach is = N n + ^h^rNr. Here ba is 
the effective branching ratio and N is the number of 
i-th hadron [49]. For details, see Refs. [46-48]. This 
particular model does not reproduce ^(piS^) nor px 
spectra at RHIC as shown in Fig. 20 in Ref. [41]. Note 
that a model is called "chemical freezeout (CFO)" when 
the number of hadrons Ni instead of N is conserved. 
This means even inelastic scatterings through strong 
interaction cease to happen. 

Hadronic cascade model (EG). One can utilize a 
hadronic cascade model just after the phase transition 
between the QGP and hadron phases [50, 51]. This ap- 
proach dynamically describes both chemical and thermal 
freezeouts without assuming explicit freezeout tempera- 
tures. Viscous effects are effectively taken into account 
through the non-zero mean free path among the particles 
(see, e.g. Eqs. (1) and (9)). 



Hydrodynamic results from the above three classes of 
hadro-dynamical models are summarized in Table I. For 
reviews of hydrodynamic results at the RHIC energies, 
see also [52-54]. As long as the space-time evolution of 
the hadron phase is concerned, the approach HC seems 
to be the most realistic one among the three classes. The 
second best model should be the model PCE from the ex- 
perimental data of particle ratios and spectra. The mod- 
els CE and HC reproduce vi(j)t) for pions and protons in 
low pt regions, whereas the model PCE fails. The fail- 
ure of PCE for vi(j>t) is particularly perplexing since the 
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TABLE I: Summary of hydrodynamic results. Hydrodynamic 
results of V2(pr) and px spectra for pions and protons are 
compiled in Fig. 20 in Ref. [41]. CE, PCE, and HC stand for, 
respectively, chemical equilibrium, partial chemical equilib- 
rium, and hadronic cascade. Thermal freezeout temperature 
is assumed to be chosen to reproduce particle spectra in the 
model CE. Note that the model CE can reproduce the shape 
of dN /prdpr for protons, not its yield. 

Observables model CE model PCE model HC 



V2{pt\ m) 


yes 


no 


yes 


dN/prdpr 


yes 


no 


yes 


ratios 


no 


yes 


yes 



spatial azimuthal asymmetry is mostly gone by the time 
hadronization is competed. In the following sections, we 
shall show why the differences between the models CE 
and PCE appears and why the result from the perfect 
fluid model CE eventually looks similar to the one from 
the highly dissipativc model HC. The key quantities to 
understand these differences are found to be the tempo- 
ral behavior of the mean transverse momentum/energy 
and the ratio of the particle number to the entropy. 



IV. TEMPORAL BEHAVIOR OF TRANSVERSE 
ENERGY 

In this section, we briefly review the time evolution of 
total transverse energy and transverse energy per particle 
in relativistic heavy ion collisions within a framework of 
the Bjorkcn solution for longitudinal expansion [22]. 

Assuming the Bjorkcn scaling solution, the time evo- 
lution of the entropy density is represented by s(r) = 
sqTq/t. Here so is the initial entropy density and To is 
the initial time. As long as a perfect fluid is considered, 
the entropy per unit rapidity is conserved 



dS 
dy 



— = I d 2 x±Ts(r) = Aj_t s 0i 



(14) 



where A± is the transverse cross section of a nucleus. 
Here we assume a smooth function for the equation of 
state (EOS). For EOS with P = c 2 s t (0 < c 2 < 1/3), the 
time evolution of energy density becomes 



e(r) 



(?) 



l+cf 



(15) 



where eo is the energy density at to. Thus the transverse 
energy per unit rapidity 



dE T 
dy 



d 2 xj_Te(r) = Aj_e Q T ( — 



(16) 



decreases with proper time due to PdV work in spite 
of the conservation of entropy [55, 56]. In the following 
discussion in this section, we consider three EOS models 
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FIG. 3: The ratio of number density and entropy density for 
massive pions. Solid (Dashed) line represents the ratio for 
chemical equilibrium (chemical freezeout). Chemical freeze- 
out temperature is assumed to be T ch = 170 MeV. 



for pions and see time evolution of total transverse energy 
and transverse energy per pion. 

Massless Pions. The number density n is propor- 
tional to the entropy density s for massless pions: n = 
(d/n 2 )((3)T 3 , s = d(47r 2 /90)T 3 . Inserting d = 3 and 
£(3) ~ 1.2, one obtains s f=a 3.6n. From Eq. (14), the 
number of pions per unit rapidity dN/dy is also con- 
served. Therefore the mean transverse mass, which is 
identical to the mean transverse momentum in the mass- 
less pion case, (tot) = (dEx / dy) / (dN / dy) oc (dEx/dy) 
decreases with t from Eq. (16). 

Massive Pions in Chemical Equilibrium. The propor- 
tionality between the number density and the entropy 
density is approximately valid only for ultra-rclativistic 
particles (T 3> to). In relativistic heavy ion collisions, 
the typical temperature is around the order of the pion 
mass. So pions are no longer relativistic particles in this 
particular situation. The number density and the entropy 
density for pions are evaluated in the usual prescription 
of statistical physics. Thus the ratio of the number N and 
the entropy S increases with temperature due to the finite 
mass of pions as shown in Fig. 3. Note that the volume 
V is canceled and that N/S equals to the ratio of their 
densities n/s. Therefore (tot) can increase with proper 
time (or with decreasing temperature of the system) even 
as dEx/dy = (mx)dN/dy decreases from Eq. (16). This 
"local reheating" can occur because the total transverse 
energy is distributed among a smaller number of pions at 
lower temperature. The resulting temporal behavior of 
(tot) is determined through competition of how rapidly 
dEx/dy and dN/dy decrease with proper time. As we 
will see in the next section, (px) increases with proper 
time in hydrodynamic simulations with chemical equilib- 
rium EOS. 

Chemically Frozen Massive Pions. When the system 
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expands strongly, inelastic collisions can cease to happen. 
So one can think about a situation in which the system 
keeps only thermal equilibrium through clastic scatter- 
ing. Analyses based on statistical models and blast wave 
models show that thermal freezeout temperature T th is 
smaller than chemical freezeout temperature T ch at AGS, 
SPS and RHIC energies. Moreover, T ch is found to be 
close to the (pseudo) critical temperature T c . This in- 
dicates that the hadron phase in relativistic heavy ion 
collisions is in thermal equilibrium, not in chemical equi- 
librium [42, 43]. Usually, the term "chemical equilib- 
rium" is associated with a state equilibrated among a 
finite number (> 1) of compositions in the system. Here 
we simply use the term "chemical freezeout" in spite of 
one hadronic species, i.e. pions. This means the num- 
ber of pions per unit rapidity is fixed below T ch . The 
entropy is also conserved as long as a perfect fluid is con- 
sidered, so the ratio of the number density and the en- 
tropy density is a constant of motion similar to the case 
for massless pions. It is interesting to mention that the 
entropy is not generated even in the evolution of chem- 
ically frozen fluids. Entropy production originates from 
the source term in the rate equation for chemical non- 
equilibrium processes. The number of hadron is, how- 
ever, conserved after chemical freezeout. It is easy to 
show the conservation of entropy d^S* 1 = from the 
conservation of energy and momentum d^T^ = 0, where 
T^ u = (e + P)u^u u - Pg^ v , and the conservation of the 
number of hadrons d^N- 1 = 0. In this sense, one needs 
to distinguish "chemical freezeout" from "chemical non- 
equilibrium" . In the chemical non-equilibrium process, 
the system is approaching to chemical equilibrium state, 
i.e. the maximum entropy state through inelastic scatter- 
ing. Entropy is certainly generated during this process. 
Contrary to this, chemical freezeout means that the sys- 
tem suddenly leaves from the chemical equilibrium state 
by keeping particle ratios due to the strong expansion. 

Figure 3 shows comparison of the ratio of pion number 
N and its entropy S between chemical equilibrium EOS 
and chemically frozen EOS. Here chemical freezeout tem- 
perature is taken as being T ch = 170 MeV. Similar to the 
massless pion case, (my) in the chemical freezeout case 
decreases with decreasing decoupling temperature. As 
long as the Bjorken scaling solution for longitudinal ex- 
pansion is considered, transverse expansion does not spoil 
the above discussion: PdV work done in the transverse 
direction is simply converted into the kinetic energy of 
fluid elements. The resultant slope of pr spectrum for 
pions should become softer at lower decoupling (thermal 
freezeout) temperature. Quantitatively, the px slope is 
insensitive to the choice of T th since dEr/dy oc t~ Cs 
changes only gradually in the late stage. The universal 
behavior of the pt slope is already confirmed in the hy- 
drodynamic simulations with chemically frozen (or par- 
tial chemical equilibrium) EOS [46, 48] and will be men- 
tioned in the next section. 

From the above consideration, the key quantity which 
governs the transverse dynamics, particularly the time 



evolution of mean transverse mass, within ideal hydro- 
dynamics is found to be the ratio of the number N and 
the entropy S. It is commonly expected that the behavior 
of the mean transverse energy/momentum is governed by 
the stiffness of the EOS. But it is not always true since the 
sound velocity of chemical freezeout EOS (or simply the 
slope of P as a function of e) is almost the same as that of 
chemical equilibrium EOS as shown in Ref. [46]. Inter- 
estingly, whether the mean transverse energy increases 
or decreases with the proper time is governed by N/S 
and the longitudinal work, not the stiffness of EOS. We 
will also mention these hydrodynamic results in the next 
section. 

To summarize this section, (mx) decreases with proper 
time for massless pions or chemically frozen massive pi- 
ons, while it can increase for massive pions in chemical 
equilibrium. We emphasize here that these conclusions 
are obtained from quite basic assumptions: the first law 
of thermodynamics (PdV work) and the Bjorken scaling 
solution. 



V. RESULTS FROM HYDRODYNAMIC 
SIMULATIONS 

We compare the hydrodynamic results from the model 
PCE with the ones from the model CE with respect to 
the EOS, space-time evolution, pr spectra, and V2{pt)- 
Hydrodynamic simulations have already been performed 
for Au+Au collisions at y / s]vw = 130 GeV and the es- 
sential results have already been obtained in Ref. [46] . In 
this section, we make an interpretation of these hydro- 
dynamic results obtained so far. In particular, we em- 
phasize that the temporal behavior of the mean p? is the 
key to understand the difference of the results between 
these two models. For further details of the hydrody- 
namic model, see also Ref. [46]. 



A. Equation of state and space-time evolution 

Chemical freezeout does not change EOS, i.e. pressure 
as a function of energy density P(e), so much [46, 47]. 
This means that the difference of chemical composition 
in the hadron phase does not affect the space-time evolu- 
tion of energy density significantly. However, at a fixed 
temperature, the energy density in chemically frozen (or 
partial chemical equilibrium) hadronic matter is consid- 
erably larger than the one for hadronic matter in chem- 
ical equilibrium. This is due to the fact that the large 
resonance population keeps in the system during the ex- 
pansion after chemical freezeout and that the mass terms 
significantly contributes to the energy density. Therefore 
the space-time evolution of temperature field does change 
significantly while P(e) remains essentially unchanged. 
The temperature of the chemically frozen system cools 
down more rapidly than that of the chemical equilibrium 
system. This leads to the reduction of life time of fluids, 
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FIG. 4: Average transverse momentum for pions as a function 
of thermal freezeout temperature for the models CE (solid) 
and PCE (dashed). 



FIG. 5: V2 for pions as a function of T th for the models CE 
(solid) and PCE (dashed). 



radial flow, and HBT radii (i?i on g and i? ou t) [46]. Lon- 
gitudinal size i?iong, which relates with the lifetime of a 
fluid through the gradient of longitudinal flow velocity 
(dvz/dz)^ 1 w Tf [57], can be interpreted by the effect of 
chemical freezeout on the life time of a fluid in hydrody- 
namics. Nevertheless, i? s idc and i? ut are still inconsistent 
with data. 



B. pt spectra and elliptic flow 

In hydrodynamic simulations with chemical equilib- 
rium, thermal freezeout temperature T th is an adjustable 
parameter. In order to fix T th , one usually utilizes the 
experimental data of pt spectra for hadrons. Reduc- 
tion of T th leads to generation of additional radial flow. 
Generically, the resulting px spectra at T th = T\ becomes 
harder than the ones at T th = Ti > T\ even though tem- 
perature, i.e. the inverse slope of the momentum distri- 
bution in the local rest frame decreases. Contrary to this 
behavior, when chemical freezeout is appropriately taken 
into account in hydrodynamic simulations, the pt slope 
becomes insensitive to T th compared to the one in the 
model CE [46]. 

To confirm these behaviors, we perform hydrodynamic 
simulations again for a particular choice of impact param- 
eter b = 5 fm and see the average transverse momentum 
(pt)- Details of the hydrodynamic models used here can 
be found in Ref. [46, 58]. Figure 4 shows (px) for pions 
as a function of T th at midrapidity y = 0. In chemi- 
cal equilibrium, (pt) increases with decreasing T th . On 
the contrary, (j>t) gradually decreases with decreasing 
T th when early chemical freezeout is taken into account. 
Even in the case that a full hydrodynamic simulation 
without boost invariant ansatz is performed and that the 
contribution from resonance decays is included, the tem- 



poral behavior of (pr) as already discussed in Sec. IV is 
seen in Fig. 4. It should be emphasized here that increase 
of (jpt) in chemical equilibrium is a direct consequence of 
neglecting chemical freezeout in hydrodynamic calcula- 
tions and, more definitely, of the experimental results of 
particle ratios. One has made full use of this unrealistic 
behavior to reproduce px spectra at the cost of hadron 
ratios in the conventional hydrodynamic calculations. 

In chemical equilibrium, the slope of elliptic flow pa- 
rameter dv2(pT)/dpT for pions is insensitive to T th and 
is consistent with the experimental data. See also Figs. 8 
and 11 in Ref. [46]. This is apparently plausible since the 
elliptic flow is a self-quenching phenomenon and is sensi- 
tive to the early stage of the collisions [59] . On the other 
hand, dv2{pr) / dpr for pions in the model PCE increases 
with decreasing T th and is ended with overprediction of 
the experimental data when T th is chosen so as to repro- 
duce the proton px spectrum and V2(pt) in the low px 
region. This means that ^(pt) for pions varies in the 
late hadronic stage (r > 10 fm/c). 

We have to be careful in understanding the difference 
between integrated elliptic flow V2 and differential elliptic 
flow V2{pt)- V2 reflects the momentum anisotropy of bulk 
matter, while V2(pt) represents how total V2 distributes 
among various particles with different px- As shown in 
Fig. 5, the integrated V2 varies only weakly with decrease 
of T th (and hence proper time t) in both cases. This is 
consistent with the self-quenching picture of elliptic flow 
again. 

The slope of i>2(pt), on the other hand, can be eval- 
uated approximately by dv2(pT)/dpT ~ v^/ipr) since 
V2(pt) for pions is almost a linear function of px [60]. In 
chemical equilibrium (CE), the moderate increase of V2 
(-13% increase as T th decreases from 160 MeV to 100 
MeV) is approximately canceled by the increase of (px) 
(~22% from 160 MeV to 100 MeV). Eventually, the ra- 
tio remains almost constant (or even decreases slightly) 
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FIG. 6: V2/(pt) for pions as a function of T th for the models 
CE (solid) and PCE (dashed). 



as shown in Fig. 6. The CE predictions for the differen- 
tial elliptic flow work because without chemical freezeout 
the slope of V2{pt) stalls accidentally and reproduces the 
experimental data. 

On the contrary, in PCE the numerator vi increases 
by ^20% and the denominator (j>t) decreases by ~10%. 
The resultant ratio i>2/(pt) therefore increases with de- 
creasing T th as shown in Fig. 6. This is the reason why 
the slope of V2(pt) in PCE increases even in the late stage 
after the spatial azimuthal asymmetry has reversed sign. 

It is now easy to understand why V2(j>t) at the SPS 
energies is so close to the one at the RHIC energy (see, 
e.g. Fig. 17 in Ref. [61]). This is due to the correlated 
change of both and (pr) from SPS to RHIC energies: 
The increase of the average V2 is compensated for by 
the increase of (pr)- The slopes of i^Cpt) therefore vary 
surprisingly weakly with the beam energy. 

Summarizing the discussion, differential elliptic flow is 
sensitive to the late hadronic stage in ideal hydrodynamic 
calculations with early chemical freezeout since the slope 
of V2{pt) is determined by the mean px, i.e. radial flow. 
The apparent consistency of V2(pr) between RHIC data 
and results based on conventional CE hydrodynamic sim- 
ulations is therefore fortuitous. 



VI. ANALYTIC APPROACH 

In order to understand the effect of chemical freezeout 
on the transverse momentum dependence of elliptic flow 
analytically, we employ the blast wave model discussed in 
Ref. [62] . In the framework of the blast wave model, one 
can choose the radial flow parameter and temperature in- 
dependently to reproduce the slope of pr spectra. How- 
ever, these two parameters are correlated in actual hydro- 
dynamic simulations. After chemical freezeout, the mean 
Pt decreases with decreasing T th , i.e. d(pr)/dT > 0, as 



already discussed in the previous sections. On the con- 
trary, the mean pt increases with decreasing T th , i.e. 
d(pT)/dT < 0, in chemical equilibrium. So we can con- 
strain the average flow velocity as a function of temper- 
ature through the condition d(pT)/dT = 0. We call the 
obtained radial flow the critical radial flow, u£ rlt (T). The 
critical radial flow ensures that the mean pt is a constant 
of motion. Qualitatively, v r (T) < uJ: rlt (T) corresponds to 
the chemically frozen system, while v r {T) > u£ nt (T) cor- 
responds to the system under chemical equilibrium. In 
the following in this section, we assume only pions which 
are dominant particles in a fluid element. 



A. Momentum Distribution 

The invariant momentum spectrum is given by the 
Cooper-Frye formula [63] in the Boltzmann approxima- 
tion: 



E 



dN 



fJ 



d 3 p (2ir) 3 



p^dun exp — 



/' 



T 



(17) 



Here g is the degree of freedom, da 11 is the element of 
freezeout hypersurfacc. p^ is the four momentum mea- 
sured in the laboratory system, p and T are, respectively, 
chemical potential and temperature at freezeout. Four 
fluid velocity (u^u^ = 1) can be parametrized as 

= cosh / o(cosh?7/, sinhpeos </>, sinhpsinc/), sinh?y^). 

(18) 

Here p is the transverse rapidity and rjf is the longitudi- 
nal rapidity which is to be identified with the space-time 
rapidity rj s in the Bjorken boost invariant solution [22]. 
One can also write cosh p = s/l + Uj_ and sinh p = u± , 
respectively. Energy of a particle in the local rest frame 
becomes 



= m T cosh(j/ - r]f)Jl + u] 



Ptu± cos{ 



Here y is the momentum rapidity and 4> p is the azimuthal 
angle of the momentum. According to Ref. [62], we also 
make the same ansatz u±(4>) = u±(l + ecos2cj)) for az- 
imuthal dependence of radial flow and take terms up to 
the first order in e 



= m T cosh(y-?7 / )y / l + 



1 



1 



-£ COS 2(j) 



— ptu± cos((p p — cj))(l + ecos2(/)). (19) 

Assuming the matter suddenly freezes out at t/, 

p^da^ = EdV = nix coshy x rdrd<\>Tjdr\ s . (20) 

Note that w_i_(</>) is supposed to include all possible az- 
imuthal anisotropic effects and that <fi dependences of T 
and p, are neglected as in the hydrodynamic approach. 
Then Eq. (17) reduces to 



dN 



mxdmxd(f> p dy 



oc 



brriT cosh ye / 



(21) 
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T 

p T u ± cos((^p - 0) ^ 

r t 

niTU 2 ^ cosh(y — rjf) cos 2</> 
Ptu± cos 20 cos(0 p — 0) 



Here. 



(22) 



B. Azimuthal Anisotropy 



By using the above momentum distribution, one can 
calculate v 2 {mT) (or V2(pt)) 
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Note that one can replace a higher order modified Bcssel 
function with lower order functions. 



Thus we obtain the same equation as Eq. (33) in Ref. [62] 



C. Incorporation of Transverse Dynamics 
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(25) 
(26) 

(27) 



Here if, and li are modified Bcssel functions. It is al- 
ways understood that the argument of Ki (Ii) is z E = 
mry/l + u]_/T (z v = pru±/T). Detailed calculations 
can be found in Appendix A. 

One can also obtain an analytic expression for the slope 
of v 2 (pt) 



From discussion in Sec. IV, we find ^(tot) /dT < 
(or d{m,T)/d(3 > 0) for chemical equilibrium pions and 
d{rriT) /dT > (or d(mr) /df3 < 0) for chemically frozen 
pions. These features are quite generic for ideal Bjorken 
fluids of pions. Obviously, the analytic approach is just 
a parametrization at frcezcout and contains almost no 
information about the time evolution of the system. For 
example, the analytic approach does not tell us anything 
about how hi increases with decrease of temperature. In 
this subsection, we try to give a dynamical meaning to 
the blast wave approach discussed in the previous sub- 
section. 

The transverse mass distribution in the analytic ap- 
proach is (see Eq. (25)) [64] 
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Thus one obtain the mean transverse mass 
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and its derivative with respect to the inverse temperature 
d(rriT) 1 
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The numerator of Eq. (37) reduces to 
[numerator of Eq. (37)] 

dmi J dm2m\m 2 l {m2 — mi)(KiIo) |, 

rriT cosh p + m<r sinh p— (3 ) KqIq 
dj3 J 



tanh p I K\I[ 



1 dp 

j3 + df3 

p T sinhp + p T coshp^/3 ) 



(40) 



The second term in the square bracket vanishes due to 
the antisymmetry (mi <-> 7712). Thus, 



numerator of Eq. (37)] 

dmi J dm 2 m\m\(rri2 - mi)Ki(zE,i)h{z p ,i) 
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The 



sign of d(mj-) /dp 1 is determined by this equation. One 
can obtain the relation between p and (3 by solving an 
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FIG. 7: Temperature dependences of radial flow in the hadron 
phase. Solid line shows a solution of the critical flow (see text 
for details) with a given initial condition (T, v r ) = (170, 0.25). 
Examples for the model CE (a = 1.2, dashed line) and CFO 
(a — 0.6, dotted line) are shown. 



equation d(rriT)/df3 = 0: 
dp , 



^45 sinh P ~ G (p) cosh P\ 
[F(p) cosh p — G(p) sinh p] = 0, 



(42) 



where, 



F(p) = J dm\dm2(m2 — m\)m\m\ 

x K 1 (z jE , 1 )/ (z p , 1 )A' (z jE , 2 )/o(2p,2), (43) 
G{p) ~ I dmidm2(rn2 — m\)m\m\\ m\ — m 2 



x K 1 {zE,i)Io{z p ,x)Kx(zE,2)h{z Pt 2). (44) 

The temperature dependence of transverse rapidity p = 
/?(/?) is obtained for a given "initial" condition ([3q,Po)- 
This particular radial flow ensures the mean transverse 
mass becomes a constant and is an upper limit of av- 
erage radial flow in chemical freezeout EOS for mas- 
sive pions. We call this solution the critical radial flow 
v"' 1 * = tanh p((3) . One can parametrize the temperature 
dependence of radial flow by introducing a parameter a 
within the analytic approach which embodies the trans- 
verse dynamics of the chemically frozen/equilibrated pion 
fluid: 



v r (T) = v r (T ch ) + a[v™(T) - v r (T ch )} 



(45) 



where v r (T ch ) = tanhp (/3o) is an initial condition for 
Eq. (42). Although the exact value of a needs much 
more involved dynamical calculation, radial flow quali- 
tatively corresponds to the chemically frozen system for 
< a < 1 and to the chemical equilibrium system for 
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result at T = T ch = 170 McV. The overall slope up 
to 1 GeV/c is gradually increasing with decreasing the 
temperature in the analytic model CFO, while the slope 
is decreasing in the analytic model CE. These results 
clearly show that the slope of V2{pt) can vary in the late 
hadronic stage and that the temporal behavior of average 
transverse momentum (pr) and radial flow v r is the key 
to understand the shape of ^(pr)- 



VII. CONCLUSION AND OUTLOOK 
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FIG. 8: Elliptic flow parameter as a function of transverse 
momentum within the analytic approach for pions. Thick 
solid line shows V2(pt)/v2 at chemical freezeout. Thin solid 
(dotted) line shows the result from the analytic model CFO 
(a = 0.6) at T th =100 (140) MeV. Dashed (dashed-dotted) 
line shows the result from the analytic model CE (a = 1.2) 
at T th =100 (140) MeV. 



a > 1. It should be mentioned here that the tempera- 
ture dependence of average radial flow can be described 
to some extent without solving hydrodynamic equations. 
Note that a should be taken as being a moderate value 
so that the total energy of the system (per unit rapidity) 
does not increase due to generation of radial flow. 

Figure 7 shows temperature dependences of radial flow. 
The solid line shows the critical radial flow i£ nt which re- 
sults from d(rriT)/dT = 0. This is obtained by solving 
Eq. (42) with an initial condition (T,v r ) = (170,0.25) 
which is consistent with a value at RHIC energies [46]. 
The critical flow distinguishes the system of chemical 
equilibrium from that of chemical freezeout: v r (T) for 
CE (CFO) should be located above (below) this line since 
d(m T ) /dT < for CE (d(m T ) /dT > for CFO). Dashed 
line (a = 1.2) shows an example of radial flow in the an- 
alytic model CE, while dotted line (a = 0.6) shows the 
one in the analytic model CFO. These results look very 
similar to the results from real hydrodynamic simulations 
as shown in Fig. 5 in Ref. [46]. 

By using these profiles for radial flow, we calculate 
V2(pt) for pions below the chemical freezeout tempera- 
ture. Hydrodynamic analysis tells us that the integrated 
t>2 is saturated within first 3-4 fm/c just after the colli- 
sion and insensitive to the late hadronic stage [53] . How- 
ever, within our analytic approach, there is no dynamical 
mechanism which saturates the integrated v-i in the late 
hadron stage. Therefore ^(pt)/^ is the quantity to be 
compared with the results from full hydrodynamic simu- 
lations. 

In Fig. 8, V2(pt)/v2 for the analytic models CE and 
CFO are represented. The thick solid line shows the 



In this paper we showed that the differential elliptic 
flow observable V2(Pt), which is a critical component 
for the interpretation of RHIC data in terms of perfect 
fluidity of the sQGP core, is sensitive to the degree of 
hadro-chemical equilibrium in late time evolution of the 
hadronic corona. If local equilibrium hydrodynamics is 
applied to the hadronic corona below T c , an inevitable 
logical impasse arises when confronting all the data on 
(1) hadron abundances (2) radial flow and (3) differen- 
tial elliptic flow. In CE hydrodynamics (2) and (3) can be 
reproduced at the expense of (1). In PCE (1) is enforced 
at the expense of (2) and (3) as summarized in Table I. 
We presented a simple analytic blast wave model to ex- 
plain these nonintuitive consequences of hadro-chemical 
(non)equilibrium in (P)CE implementations of hydrody- 
namics. 

In CE both the average transverse momentum per 
hadron (j>t) and average i>2 increase with proper time 
in the hadronic phase in a way that accidentally pre- 
serves the slope of differential elliptic flow dv2 {pt)I dpr ~ 
V2/(pt) in agreement with the data. In PCE, (px) de- 
creases due to the basic Bjorken longitudinal cooling. 
The main point is that in PCE the hadronic yields are 
fixed at T ch and the compensating CE "local reheating" 
mechanism (the conversion of heavy resonance mass back 
into internal energy which "mimics" a sort of dissipative 
effect) is absent. This is why PCE fails to describe the 
proton radial flow data. In addition, the slight increase 
of the average V2, as in CE, with proper time cannot be 
compensated for in PCE. Therefore, the slope of differ- 
ential elliptic flow dv2{pT)/dpT ~ V2/(pt) continues to 
grow in PCE during the hadronic phase, which leads to 
disagreement with RHIC data. 

The subtle interplay among (1) longitudinal expansion 
work, (2) maintenance of hadronic abundance yields, (3) 
the long time development of radial flow, and (4) the 
differential azimuthal asymmetric elliptic flow provides 
a formidable dynamical constraint on the dynamics of 
the hadronic corona. Only by abandoning ideal hydro- 
dynamics in the hadronic corona, have nonequilibrium 
hadron cascade (HC) models been able to deal with the 
interplay of the above hadron dynamics in a way consis- 
tent with present RHIC data. As discussed in Sec. II, this 
approach is natural since the viscosity to entropy ratio in 
a hadronic resonance gas below T c is too large to support 
even local thermal equilibrium. By relaxing both thermal 
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and hadro-chcmical equilibrium assumptions, the hybrid 
QGP hydrodynamics plus hadron cascade model in [51] 
has been able to account for all three major low pr ob- 
servables as summarized in Table I. The effect of viscosity 
in the hadron phase [23] substitutes for the "local reheat- 
ing" in the CE model and compensates the small growth 
of the average V2 in PCE to preserve the slope of V2(pr) 
for pions. The slope of ^(pt) is thus found to stall at 
the SPS energy from the hybrid model analysis [47]. In 
the classical transport approach, both (pr) [65] and V2 
[66] do not vary significantly when the mean free path 
among the particles becomes comparable with the typ- 
ical gradient length scales. Moreover, the shear viscous 
effect changes the momentum distribution function [67] 
and reduces the slope of V2{pr) slightly [23]. These are 
the reasons why the slope of V2(pt) does not changed 
significantly in the hadronic transport stage. 

So how robust is the statement that hydrodynamic de- 
scription at RHIC works remarkably well? We empha- 
sized that the behavior of V2 differs from that of V2{pt)- 
The integrated elliptic flow does not develop so much 
in the late hadronic stage in which cither the inviscid, 
chemical (non-) equilibrium fluid or the dissipative gas is 
assumed, whereas the differential elliptic flow depends 
largely on these assumptions. The large magnitude of 
integrated V2 observed at RHIC is reproduced only when 
a small rj/s is assumed [18]. Therefore the large V2 devel- 
oped in the early stage is obtained from the evolution of 
the sQGP core, which as discussed in Sec. II must have 
near minimal viscosity ?7sym ~ T 3 . Even though the 
minimal sQGP viscosity is larger than the viscosity of the 
HRG corona, the core exhibits near perfect fluid behavior 
due to the deconfmcment of almost all the QCD degrees 
of freedom. The near perfect fluidity of the sQGP core is 
therefore a signal of dcconfincmcnt. On the other hand, 
the breaking of local and hadro-chemical equilibrium in 
the hadronic corona is critical for this interpretation of 
RHIC data. If inviscid ideal hydrodynamics were valid in 
both sQGP and HRG phases, the crucial ^(pt) would 
be sensitive to the hadronic thermal freezeout dynamics 
and not only to the equation of state of sQGP matter. 

Perhaps most surprising in connection with the hy- 
drodynamics robustness question is the important role 
played by hadro-chemical freezeout at T ch sa T c that is 
implied by the extensive systematics of observed hadron 
abundances [44]. Without this constraint the different 
hadro-chemical results with CE and PCE would preclude 
a conclusion about the perfect fluidity of the sQGP core 
as well as the highly dissipative, imperfect fluidity of its 
hadronic resonance corona. 

Despite the success of the hybrid HC approach, there 
exist open technical questions that must be still inves- 
tigated. One important issue is the violation of energy- 
momentum conservations at the boundary between the 
QGP and hadron phases [68]. The Cooper- Frye prescrip- 
tion [63] is employed to obtain the particle distribution 
just after hadronization which is to be used as an initial 
condition in the sequential cascade calculation. The vio- 



lation is expected to be small when radial flow is large. 
Nevertheless there always exists in the space-like hyper- 
surface da^ in-coming particles which contributes to the 
number of particles negatively. This negative contribu- 
tion is omitted in the actual calculations, which causes 
the violation of energy-momentum conservation. Proper 
treatment of the boundary condition may lead to change 
the dynamics in the QGP phase [68]. In this connec- 
tion, the approximate continuity of the viscosity from the 
sQGP to the HRG phase discussed in Sec. II minimizes 
this interface problem since there is no discontinuity of 
the stress tensor including the viscous correction at T c . 

Another important future problem is the rapidity de- 
pendence of elliptic flow. The (3+l)-dimensional ideal 
hydrodynamic calculations [8, 46] have not been able 
to reproduce the observed pseudorapidity dependence 
of V2 in forward rapidity region at RHIC [3]. The for- 
ward rapidity region at RHIC is similar to the midra- 
pidity region at SPS in the sense that local particle den- 
sity (l/S)dN c \ l /dy is similar. Heinz and Kolb [69] pro- 
posed a "thermalization coefficient" to describe enhanced 
nonequilibrium effects in the low particle density region 
defined from the experimental data 1)2/ £ as a function 
of (1 / S)dN c h/ dy [4]. To address correctly the rapidity 
and beam energy dependence taking into account the 
highly viscous nature of the hadronic corona, a new hy- 
brid model must be developed in which the (3+l)-D hy- 
drodynamic model for the sQGP core is combined with 
the HRG transport approach for the dissipative hadronic 
corona. Hirano and Nara [32, 58, 70] have already devel- 
oped a dynamical model to describe three different as- 
pects of relativistic heavy ion collisions in one consistent 
framework: Color glass condensate initial conditions for 
high energy colliding nuclei, hydrodynamic evolution in 
3D space for long wave length components of produced 
matter, and quenching jets for short wave length com- 
ponents. A further unified framework by combining a 
hadronic cascade model with the above one [71] will be 
necessary to understand more quantitatively the dynam- 
ics and the properties of QCD matter produced in rela- 
tivistic heavy ion collisions at all SPS, RHIC and LHC 
energies. 

We emphasize that continued work toward such a uni- 
fied dynamical framework will be essential to further 
test our physical interpretation of RHIC data - as out- 
lined in the introduction and analyzed in sections II- VI 
- that "perfect fluidity" of the higher viscosity sQGP 
core and "imperfect fluidity" of the lower viscosity HRG 
corona taken together with hadro-chcmical equilibrium 
near T c ~ 160 — 170 MeV already provide a compelling 
set of signatures for QCD deconfmcment at RHIC. 
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APPENDIX A: DERIVATION OF EQ. (24) 

Let us recall the formulae for modified Bcssel functions 

/oo 
dxcosha;exp(— zeoshec), (Al) 
-oo 

r 

271-/2 (z) = / d</>cos20exp(zcos0). (A2) 

J — 7T 

The numerator of Eq. (23) becomes 

d<p p cos 2(f) p J d(f> J dy cosh y exp( A) 
d4>1K x (B) x 2n cos 2<f)I 2 (C) , (A3) 

where, 
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We here assume e is small, expand the modified Bessel 
function, and take the first order term with respect to s, 
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Let us also recall some useful formulae for the derivatives 
of modified Bcssel functions, 



-#„__(«) - -K n (z) 



= ^[In-l(z)+I n+1 (z)]. 

Then the numerator of Eq. (23) is proportional to 
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Here the arguments of K and / are, respectively, zg = 
iriT yl + u^/T and z p = ptu±/T. An analogous calcu- 
lation can be done for the denominator of Eq. (23). The 
result is proportional to 4-7T 2 x 2K\Iq. 
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